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B13 is a ceramide analogue and apoptosis inducer with potent cytotoxic activity. A series of 
arylpropyl sulfonamide analogues of B13 were evaluated for their cytotoxicity using MTT assays in 
prostate cancer PC-3 and leukemia HL-60 cell lines. Some compounds (4, 9, 13, 14, 15, and 20) showed 
stronger activities than B13 in both tumor cell lines, and compound (15) gave the most potent activity 
with IC50 values of 29.2 and 20.7 /i M, for PC-3and HL-60 cells, respectively. Three-dimensional 
quantitative structure-activity relationship (3D-QSAR) analysis was performed to build highly reliable 
and predictive CoMSIA models with cross-validated q 2 values of 0.816 and 0.702, respectively. Our 
results suggest that long alkyl chains and a 1R, 2R configuration of the propyl group are important 
for the cytotoxic activities of arylpropyl sulfonamides. Moreover, the introduction of small hydrophobic 
groups in the phenyl ring and sulfonamide group could increase biological activity. 
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INTRODUCTION 

Apoptosis is an endogenous program of cell death that 
is mediated through biochemical modulators, and it is in- 
duced by anticancer agents, tumor necrosis factor, and ion- 
izing radiation [1,2]. Because ceramide is known to be a 
key modulator of apoptosis and its accumulation in cells 
leads to apoptosis, many approaches have been explored to 
increase endogenous ceramide [3-5]. These approaches in- 
clude the application of short-chain ceramides (C6-ceramide) 
[6-8] and induction of ceramide by modulation of ceram- 
ide-metabolizing enzymes (LCL 15 and LCL 16) [9-13]. 

In recent years, ceramide analogues have emerged as a 
new strategy for cancer therapy [14-17]. The ceramide ana- 
logues, such as sphingosine, sphinganine, and C2-eramide, 
were known to induce apoptotic cell death in the malignant 
melanoma, colon cancer cells, and prostate cancer cells 
[18-21]. It was reported that the activity of ceramide ana- 
logues is influenced by the stereochemistry and chain 
length of the amide moiety. Analogues with a ceramide con- 
figuration (2S, 3R) are not active, whereas their stereo- 
isomers inhibit cell growth. The ceramides with medium 
size alkyl chains (Cn~Ci5) are active compounds, but ce- 
ramides with shorter (Ci~Cg) and longer (Ci7~C23) alkyl 
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chains are less active [22]. The aromatic analogues of ce- 
ramide like D-e-Mapp with a IS, 2R configuration increase 
endogenous ceramide to induce apoptosis in many cancer 
cell lines [23-26]. Another active analogue, B13 with a 1R, 
2R configuration, exhibited cytotoxicity in malignant mela- 
noma, colon, and prostate cancer cells [27-29]. B13 also in- 
duced apoptosis in cancer cells, increased the cellular ce- 
ramide level, inhibited metastasis, and had no effect on nor- 
mal cells [29-31]. 

It has been reported that isosteric replacement of the 
amide group in a ceramide by urea or amine can increase 
the ceramidase inhibitory activity and cytotoxicity [32]. 
And the carboxyl moiety of ceramide can also be bioisosteri- 
cally replaced by a sulfone group. Sulfonamides are known 
to have various biological activities, including hypoglycemic, 
diuretic, antibacterial activity, and cytotoxicity [33-36]. 
Extending these findings to the ceramide, we suspected 
that introduction of a phenyl ring, as in D-e-Mapp and B13, 
altering alkyl chain lengths and stereochemistry, and bio- 
isosterical replacement of the carboxyl group with a sulfone 
group would improve and modify the biological properties 
of ceramide. 

As a continuation of our previously reported work, in 
which arylpropyl thiourea analogues (Fig. 1) were studied 
[37], we determined the cytotoxic activity of 21 new ar- 
ylpropyl sulfonamide analogues (Fig. 1) in prostate cancer 
and leukemia cell lines to evaluate the effects of R and S 
enantiomers in propyl moiety, a p-nitro group in the phenyl 



ABBREVIATIONS: DMSO, dimethyl sulfoxide; D-e-Mapp, (lS,2R)-2- 
N-(tetradecanoylamino)-l-phenyl-l-propanol; B13, (lR,2R)-2-N-(tet- 
radecanoylamino) - 1 - (4' -nitrophenyl) - 1 , 3-propanediol; PBS, phosphate- 
buffered saline; MTT, 3-[4,5-dimethylthiazol-2-yl]-2,5-diphenyltetra- 
zolium bromide; 3D-QSAR, three dimensional quantitative structure 
activity relationship; CoMSIA, comparative molecular similarity 
indices analysis. 



237 



238 



YJ Hwang, et al 





Arylpropyl thioureas 




0"^-0 
Arylpropyl suffonamides 



Fig. 1. Ceramide, Mapp, B13, and 
arylpropyl compounds. 



ring, and the alkyl chain length of the sulfonamide moiety. 
We also perform QSAR analysis to investigate the relation- 
ship between the structural features and cytotoxicity of ar- 
ylpropyl sulfonamides and to develop more potent anti- 
cancer agents. 



METHODS 



Materials 



A series of arylpropyl compounds were previously syn- 
thesized in our lab. Phosphate-buffered saline (PBS) was 
purchased from Boehringer Mannheim. Dimethyl sulfoxide 
(DMSO), 3-[4,5-dimethylthiazol-2-yl]-2,5-diphenyltetrazolium 
bromide (MTT), and other reagents were obtained from 
Sigma. 

In vitro cytotoxic assay 



the CoMSIA studies. 

Molecular modeling and alignment 

The modeling software Sybyl-X 1.3 was used for the 
structure building, molecules modeling, partial least 
squares, and conducting CoMSIA [40]. The structures of 
compounds were generated with a sketch tool and energy 
minimization was performed using a TRIPOS force field 
with the Powell method and conjugate gradient termination. 
The atomic charges of molecules were calculated using the 
Gasteiger- Huckel charges. Simulated annealing was used 
to determine the low-energy conformations. One of the most 
important requirements for CoMSIA models is that the 3D 
structures of molecules are aligned to a suitable conforma- 
tional template. The molecular alignment was achieved by 
the fitting atoms method. In Table 1, the bold line repre- 
sents a common substructure and B13 was used as a tem- 
plate molecule in the alignment. 



The cytotoxicity of arylpropyl compounds was evaluated 
in two human tumor cell lines: prostate cancer PC-3 cells 
and leukemic cancer HL-60 cells. The cytotoxicity was de- 
termined using a MTT-based colorimetric assay [38]. The 
cells were treated as described in the Table 1 legend and 
the results from the assay are shown in Table 1 [39]. 

Data sets 

Twenty-one arylpropyl compounds with cytotoxic IC50 
values ranging from 20.7 to 267.3 ^M were used to carry 
out 3D-QSAR analysis. Their molecular structures are il- 
lustrated in Table 1. The test set was selected based on 
their various configurations at Ci and C2 and used for ex- 
ternal validation of the 3D-QSAR models. The training set 
was consisted of seventeen compounds including B13 and 
the test set was made of the following four different config- 
uration structures, compound 2(1R,2S), 8(1S,2R), 12(1R,2R), 
and 17(1S,2S). All IC50 values were transformed into pICso 
( — log ICbo) values and used as the dependent variables in 



CoMSIA 3D-QSAR models 

CoMSIA is based on the relationship between the bio- 
logical activity and structural properties of compounds 
when the receptor structure is not known. It evaluates the 
five physicochemical properties: steric, electrostatic, hydro- 
phobic, hydrogen bond acceptor, and hydrogen bond donor 
fields. The CoMSIA method involves a common probe atom 
and similarity indices determined in regularly spaced grid 
points for the aligned molecules. The common probe atom 
with a radius of 1.0 A, charge of +1, hydrophobicity of +1, 
hydrogen bond donating of + 1, and hydrogen bond accept- 
ing of + 1 was used to calculate the five fields. A default 
value of 0.3 was used for the attenuation factor. 

Partial least squares (PLS) analysis 

In the PLS analysis, three CoMSIA descriptors (electro- 
static, hydrophobic, and hydrogen bond acceptor fields) 
were used as independent variables, and the pICso values 
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Table 1. Structures and cytotoxic activities of arylpropyl sulfonamides 




Cytotoxicity IC50 (fM) 

Compounds R, R 2 R 3 Configuration Prosta te cancer Leukemia 

(PC-3) (HL-60) 



1 


H 


H 


C7H15 


1R, 


2S 


189.3 


129.8 


2 


H 


H 


CsH 17 


1R 


2S 


95.3 


77.6 


3 


H 


H 


C9H19 


1R, 


2S 


81.7 


53.7 


4 


H 


H 


C11H23 


1R, 


2S 


44.9 


28.5 


5 


H 


H 


C13H27 


1R, 


2S 


45.1 


34.3 


6 


H 


H 


C7H15 


IS, 


2R 


188.5 


130.4 


7 


H 


H 


CsH 17 


IS, 


2R 


80.8 


62.8 


8 


H 


H 


C9H19 


IS, 


2R 


64.3 


44.9 


9 


H 


H 


C11H23 


IS, 


2R 


40.5 


33.1 


10 


H 


H 


C13H27 


IS, 


2R 


67.8 


44.4 


11 


OH 


N0 2 


C7H15 


1R, 


2R 


98.4 


116.6 


12 


OH 


N0 2 


CsH 17 


1R, 


2R 


52.8 


54.9 


13 


OH 


N0 2 


C9H19 


1R, 


2R 


31.8 


27.6 


14 


OH 


N0 2 


C11H23 


1R, 


2R 


39.1 


24.7 


15 


OH 


N0 2 


C13H27 


1R, 


2R 


29.2 


20.7 


16 


OH 


N0 2 


C7H15 


IS, 


2S 


> 267.3 


160.6 


17 


OH 


NO2 


CsH 17 


IS, 


2S 


75.1 


67.6 


18 


OH 


N0 2 


C9H19 


IS, 


2S 


51.0 


41.4 


19 


OH 


N0 2 


C11H23 


IS, 


2S 


56.0 


37.7 


20 


OH 


NO2 


C13H27 


IS, 


2S 


35.9 


23.0 


B13 








1R, 


2R 


79.3 


33.6 



The cells were plated at a density of approximately lxlO 4 cells/well in 96-well plates. Each well contained 180 u 1 of medium and 
20/^1 of lOxconcentration of prepared compounds or PBS were added. After 96 h of culture, 0.1 mg of MTT was added to each 
well and incubated at 37°C for 4 h. The plates were centrifuged at 450xg to precipitate the formazan crystals. The medium was 
rem9oved and 150 ;jL of DMSO was added to each well to dissolve the formazan. In this assay, MTT was converted to blue formazan 
by mitochondrial dehydrogenase. The intensity of the blue color was measured with a microplate reader at a wavelength of 540 
nm. The measured mean values were expressed as the IC50, the concentration that reduced the optical density of the treated cells 
by 50% with respect to the untreated controls. 



Table 2. PLS analysis of CoMSIA 3D-QSAR models 



PLS Statistics 


Prostate 
cancer (PC-3) 


Leukemia 
(HL-60) 


q 2a) 


0.816 


0.702 


j» 


0.999 


0.997 


SEE C) 


0.009 


0.018 


F d) 


3841.601 


1000.475 


PLS component"' 


4 


4 


r pred 


0.776 


0.974 


Field Steric 






contribution Electrostatic 


0.416 


0.518 


Hydrophobic 


0.335 


0.254 


H-bond donor 






H-bond acceptor 


0.249 


0.228 



d) q 2 , cross-validated correlation coefficient from leave-one- out 
(LOO); b) r 2 , non-cross-validated correlation coefficient; c) SEE, 
standard error of estimate; *F, F-test value; e) PLS component, opti- 
mum number of components; V pre d, predicted correlation coefficient. 



were used as dependent variables. The predictive pICso val- 
ues of the models were evaluated by leave-one-out (LOO) 
cross-validation. In the LOO method, one compound is re- 



moved from the data set and its biological activity is pre- 
dicted with the model derived from the rest of data set. 
The LOO method determines the optimum number of com- 
ponents, which are then used for the non-cross-validated 
analysis. To test the utility of the model as a predictive 
tool, the test set compounds that were not used in the model 
generation were predicted. 



RESULTS 

In vitro cytotoxic activity 

The cytotoxicities of 21 compounds were evaluated in 
vitro and presented in Table 1. B13 gave moderate cytotox- 
icity with IC50 values of 79.3 and 33.6 fiM for prostate can- 
cer PC-3 and leukemia HL-60 cells, respectively. The IC50 
values of the other structures ranged from 29.2 to 267.3 u M 
for PC-3 cells and 20.7 to 160.6 for HL-60 cells. 

For prostate cancer PC-3 cells, the long alkyl chain 
(C13H27 and C11H23) compounds (4, 5, 9, 10, 14, 15, 19, and 
20) exhibited more potent activities than B13 to give IC50 
values of 44.9, 45.1, 40.5, 67.8, 39.1, 29.2, 56, and 35.9 uM, 
respectively. However the short alkyl chain (C7H15 and 
CsH 17) compounds (1, 2, 6, 7, 11, and 16) showed less potent 
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Table 3. Residuals of the predicted cytotoxicities (pIC50) of training set 



Prostate cancer cells (PC-3 cells) 



Leukemia cells (HL-60 cells) 



vumpuunua 


Experimental*' 


Predicted 1 "' 


Residual 1 ' 


Experimental a ' 


Predicted b ' 


Residual'' 


1 


q no 
a. iZ 


Q 1 A 


n no 
U.Uz 




Q OS 
d.OD 


n no 
U.IM 


Q 
O 


a no 
4.U9 


A flC 


U.U1 


4.27 


A on 
4.Z / 


U.UU 


4 


4.35 


4.34 


0.01 


4.55 


4.55 


0.00 


r 
O 


A QC 


A QC 

4.ob 


U.UU 


4.46 


4.45 


U.UZ 


6 


3.72 


3.73 


-0.01 


3.88 


3.93 


-0.04 


7 


4.09 


4.08 


0.01 


4.20 


4.20 


0.00 


9 


4.39 


4.39 


0.00 


4.48 


4.48 


0.00 


10 


4.17 


4.17 


0.00 


4.35 


4.36 


-0.01 


11 


4.01 


4.01 


0.00 


3.93 


3.93 


0.00 


13 


4.50 


4.50 


0.00 


4.56 


4.55 


0.01 


14 


4.41 


4.41 


0.00 


4.61 


4.62 


-0.01 


15 


4.53 


4.54 


-0.01 


4.68 


4.67 


0.01 


16 


3.57 


3.56 


0.01 


3.79 


3.79 


0.00 


18 


4.29 


4.28 


0.01 


4.38 


4.38 


0.00 


19 


4.25 


4.25 


0.00 


4.42 


4.43 


0.00 


20 


4.44 


4.45 


-0.01 


4.64 


4.64 


0.00 


B13 


4.10 


4.10 


0.00 


4.47 


4.48 


-0.01 


Average 






0.01 






0.01 



"'Experimental cytotoxic activity; b 'predicted activity by the CoMSIA model with electrostatic, hydrophobic, and hydrogen bond accept- 
or fields; "'difference between the experimental and predicted activities; The pICso (—log ICeo) values were converted from IC50 values. 



Table 4. Residuals of the predicted cytotoxicities (pICso) of test set 



Compounds 


Prostate 


cancer cells (PC-3 cells) 


Leukemia cells (HL-60 cells) 




Experimental'*' 


Predicted b) 


Residual"' 


Experimental"' 


Predicted b> Residual'' 


2 


4.02 


4.02 


0.00 


4.11 


4.12 


0.01 


8 


4.19 


4.17 


0.02 


4.35 


4.34 


0.01 


12 


4.28 


4.27 


0.01 


4.26 


4.23 


0.03 


17 


4.12 


4.21 


-0.08 


4.17 


4.14 


0.03 


Average 






0.03 






0.02 



Experimental cytotoxic activity; predicted activity by the CoMSIA model with electrostatic, hydrophobic, and hydrogen bond accept- 
or fields; "'difference between the experimental and predicted activities; The pICso (—log IC50) values were converted from IC50 values. 



activities. The compounds with a 1R, 2R configuration (11, 
12, 13, 14, and 15) were more active than corresponding 
compounds with other configurations. Compounds (11-20) 
with C3-OH and a para nitro group in the phenyl ring had 
similar activities to the compounds (1-10) without these 
groups except for the compounds with a 1R, 2R configuration. 

The cytotoxic activity in leukemia HL-60 cells was sim- 
ilar to that observed for prostate cancer PC-3 cells. The long 
alkyl chain compounds showed more potent cytotoxicity 
than the short alkyl chain compounds. Some of long alkyl 
chain (C11H23 and C13H27) analogues (4, 9, 14, 15, and 20) 
were more active than B13 with IC50 values of 28.5, 33.1, 
24.7, 20.7, and 23 uM, respectively. The compounds with 
a 1R, 2R configuration (11, 12, 13, 14, and 15) produced 
more cytotoxicity than corresponding compounds with other 
configurations. Compounds (4, 9, 13, 14, 15, and 20) with 
a long alkyl chain showed more potent activities than B13 
against both tumor cell lines. 



3D-QSAR analysis 

The statistical data for the CoMSIA model are summar- 
ized in Table 2. The residuals between the experimental 
and predicted cytotoxicity (pICso) are listed in Table 3 and 
4, respectively. 

Prostate cancer PC-3 cells 

The PLS analysis of the CoMSIA model with electro- 
static, hydrophobic, and hydrogen bond acceptor fields pro- 
vided a cross-validated coefficient q 2 value of 0.816 and the 
non-cross-validated coefficient r 2 value of 0.999 with the op- 
timum components number of 4, SEE of 0.009, and F value 
of 3841.601. The predictive ability of the model is expressed 
by the r 2 prel i value, which is analogous to the cross-validated 
q 2 , and the predicted r 2 value of this model was 0.776. The 
experimental pICso values, predicted pICso values, and the 
residuals for training and test sets are shown in Table 3 
and 4, respectively. Fig. 2A describes the correlation be- 
tween the experimental versus predicted pICso values for 
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4.4 - 


a 




cted 


4.1 - 


■u 

aj 






3.8 - 



3 5 



y = 0.998X + O.C 
= 0.999 




3.8 



4 : 

Experimental pICso 



4.7 



(A) 



y = 0.995X + 0.023 
R2 = 0.996 



(B) 



Experimental pICso 

Fig. 2. Graph of the correlation between experimental and predicted 
activities for training and test set compounds. (A) Prostate cancer 
PC-3 cells. (B) Leukemia HL-60 cells. The ICso values were 
transformed into pICso ( — log ICso) values (■: training set compounds, 
A: test set compounds). 



the training and test set compounds. 

Leukemia HL-60 cells 

The CoMSIA model with electrostatic, hydrophobic, and 
hydrogen bond acceptor fields gave the cross-validated q 2 
value of 0.702, non-cross-validated r value of 0.997, and 
PLS components of 4, respectively. With this model, the 
predicted r 2 value (r 2 pre d) was 0.974 and the predicted pICso 
values and residuals of the training set and test set are 
presented in Table 3 and 4, respectively. The graph of ex- 
perimental versus predicted pICso values for the training 
and test sets is shown in Fig. 2B. 



DISCUSSION 

In our in vitro biological assay, compound (15) showed 
strong activity against both prostate cancer PC-3 and leu- 
kemia HL-60 cells with IC 5 o values of 29.2 and 20.7 fiM, 
respectively. The compounds with a C13H27 alkyl chain (5, 
15, and 20) gave more potent cytotoxicity than B13, which 
also had a C13H27 alkyl chain, suggesting that the amide 
group of B13 could be replaced by a sulfonamide group to 
increase activity. This finding provides some clues to guide 
the design of new B13 analogues, since the sulfonamides 
are generally more water soluble and stable to hydrolysis 
by amidase than the amide compounds. The stereochemistry 
in the propyl group was not critical to biological activity, 
although the compounds with a 1R, 2R configuration (11, 
12, 13, 14, and 15) showed more activity than corresponding 
compounds with other configurations. In general, com- 



Fig. 3. 3D-contour maps of the CoMSIA models with electrostatic, 
hydrophobic, and hydrogen bond acceptor fields. (A) Prostate cancer 
PC-3 cells. (B) Leukemia HL-60 cells. Compound (15) is shown 
within the fields (blue, favorable electrostatic; red, unfavorable 
electrostatic; yellow, favorable hydrophobic; white, unfavorable 
hydrophobic; magenta, favorable hydrogen bond acceptor; cyan, 
unfavorable hydrogen bond acceptor). 



pounds with long alkyl chains produced increased activity 
compared to compounds with short alkyl chains. 

For prostate cancer PC-3cells, the statistical data from 
the CoMSIA model showed a high cross-validated q 2 value 
(0.816) and non-cross-validated coefficient r 2 value (0.999). 
Its q 2 value, higher than 0.5, suggested that this model was 
predictive for validation. The corresponding field con- 
tributions of electrostatic (E), hydrophobic (H), and hydro- 
gen-bond acceptor (A) were 41.6%, 33.5%, and 24.9%, 
respectively. In the CoMSIA model of leukemia HL-60 cells, 
good predictivity was represented by high q" (0.702) and 
r 2 (0.997) coefficient values. The electrostatic (51.8%), hy- 
drophobic (25.4%), and hydrogen-bond acceptor field (22.8%) 
contributions indicate that electrostatic ligand-receptor in- 
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teractions are mainly involved in the cytotoxicity. The small 
average residual value (0.01) between experimental and 
predicted activities in Table 3 shows that the predicted ac- 
tivities from CoMSIA models are correlated well with the 
experimental activities. To validate the predictivity of 
CoMSIA models, the diverse configurations at Ci and C2 
were selected for the test set and their small average re- 
sidual values (0.03 and 0.02) indicate that both models are 
accurate to predict the cytotoxicity of the test set molecules 
(Table 4). 

The CoMSIA 3D-contour maps were used to visualize the 
information of the derived 3D-QSAR models for PC-3 and 
HL-60 cells (Fig. 3). The maps use the characteristics of 
the compounds, which are important for biological activity 
and show the regions around molecules where increased or 
decreased activities are expected based on physicochemical 
property changes in the molecules. In an electrostatic con- 
tour map, the blue and red regions are favorable for positive 
and negative charges, respectively. The yellow contours in- 
dicate where hydrophobic groups (e.g. methyl and phenyl) 
increase biological activity, whereas the white contours 
show regions where hydrophilic groups (e.g., hydroxyl, car- 
bonyl, and amino) increase activity. The magenta area fa- 
vors a hydrogen bond acceptor (electronegative atoms such 
as fluorine, oxygen, and nitrogen), while the cyan repre- 
sents an unfavorable region. The molecule in the contour 
maps is compound (15), which shows the most potent 
cytotoxicity. 

The CoMSIA contour map of prostate cancer PC-3 cells 
(Fig. 3A) was similar to that of leukemia HL-60 cells (Fig. 
3B). The contour map of PC-3 cells (Fig. 3A) revealed a blue 
contour at the NH and red contours near the oxygen atoms 
of the sulfonamide, which helps to explain the strong cyto- 
toxic activity of compound (15). The yellow contours were 
shown at the ortho and meta positions of the phenyl ring, 
suggesting that the introduction of a naphthyl ring instead 
of phenyl ring or small alkyl group like CH 3 at these posi- 
tions would increase activity. This finding was supported 
by the compounds (6, 11, and 16) in which yellow contours 
were superimposed with NO2 or sulfonamide groups to give 
less potent activities. The yellow contours at the phenyl 
ring in compounds (4 and 9) and the white contours at NH 
moiety in compounds (15 and 20) provide clues about their 
high cytotoxicity. The magenta contour is near the oxygen 
atoms of the sulfonamide and the cyan contours were 
around the hydrogen atoms of the NH. These magenta and 
cyan contours appear in similar areas in compounds (4, 9, 
13, and 20), thus explaining their high potency. 

The contour map of HL-60 cells (Fig. 3B) displayed the 
blue contour at NH and red contours around the oxygen 
atoms of the sulfonamide, demonstrating the potent activity 
of the compound (15). The yellow contours were placed at 
the meta positions of the phenyl ring and C2-H. The yellow 
contours appeared around the NO2 or sulfonamide groups 
in compounds (6, 11, and 16) and the phenyl ring and meth- 
yl group in compounds (4 and 9) to explain low potency 
of the compounds (6, 11, and 16) and high potency of the 
compounds (4 and 9). The magenta contour was around the 
oxygen atoms of Ci-OH and the sulfonamide and cyan con- 
tours were near the hydrogen atoms of NH. In compounds 
(4, 9, 13, and 20), a similar distribution of magenta and 
cyan contours was found indicating their strong cytotoxicity. 

In conclusion, compounds (4, 9, 13, 14, 15, and 20) gave 
more potent cytotoxicity than B13 in both tumor cell lines. 
The good predictivity of the CoMSIA models was explained 



by high q 2 (0.816 and 0.702) and r 2 (0.999 and 0.997) co- 
efficient values. The low average residual values (training 
set for 0.01; test set for 0.03 and 0.02) in Table 3 and 4 
indicate that the predicted activities from the CoMSIA 
models of PC-3 and HL-60 cell lines correlates well with 
the experimental activities. The cytotoxic assay and 
3D-QSAR analysis suggested that the sulfonamide group, 
long alkyl chains, and 1R, 2R configuration are important 
for the cytotoxic activities. Moreover, our results suggest 
that the introduction of hydrophobic groups in the phenyl 
ring will increase the cytotoxicity. 
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